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Abstract 

Numerous non-standard dynamics are described by contact-like effective interactions that 
can manifest themselves only through deviations of the cross sections from the Standard 
Model predictions. If one such deviation were observed, it should be important to definitely 
identify, to a given confidence level, the actual source among the possible non-standard 
interactions that in principle can explain it. We here estimate the “identification” reach 
on different New Physics effective interactions obtainable from angular distributions of 
lepton pair production processes at the planned International Linear Collider with polarized 
beams. The models for which we discuss the range in the relevant high mass scales where 
they can be “identified” as sources of corrections from the Standard Model predictions, 
are the interactions based on gravity in large and in TeV -1 extra dimensions and the 
compositeness-inspired four-fermion contact interactions. The availability of both beams 
polarized in many cases plays an essential role in enhancing the identification sensitivity. 


1 Introduction 


Numerous new physics (NP) scenarios, candidates as solutions of Standard Model (SM) 
conceptual problems, are characterized by novel interactions mediated by exchanges of 
very heavy states with mass scales significantly greater than the electroweak scale. In 
many cases, theoretical considerations as well as current experimental constraints indicate 
that the new objects may be too heavy to be directly produced even at the highest energies 
foreseen at future colliders, such as the CERN Large Hadron Collider (LHC) and the e + —e~ 
International Linear Collider (ILC). In this situation the new, non-standard, interactions 
could only be revealed by indirect, virtual, effects manifesting themselves as deviations of 
measured cross sections from the SM predictions. 

At the available “low” energies provided by the accelerators, where we study reactions 
among the familiar SM particles, effective “contact” interaction Lagrangians represent the 
convenient theoretical tool to physically parameterize the effects of the above mentioned 
non-standard interactions and, in particular, to test the corresponding virtual high mass 
exchanges. Clearly, in this framework the transition amplitudes are parameterized as power 
expansions in the (small) ratios between the Mandelstam variables of the process under 
study and the high mass scales squared. The sensitivity to the searched for signals will 
therefore be increased by the colliders high energy (and high luminosity). 

Since, in principle, the observation of a correction to the SM cross section may by itself 
not enable to unambiguously identify its source among the different possible explanations, 
suitable observables with enhanced sensitivity to the individual non-standard scenarios 
and/or convenient statistical criteria in the data analysis must be defined in order to dis¬ 
criminate the nature of the relevant new physics model against the others. As regards the 
search of effective contact interactions at the planned high energy colliders, it should there¬ 
fore be desirable to assess for each non-standard scenario not only the “discovery reach”, 
i.e., the maximum value of the relevant mass scale below which it produces observable 
corrections to the SM predictions, but also the upper limit of the range of mass scale val¬ 
ues where the considered scenario not only produces observed deviations, but can also be 
discriminated from the other potential sources of the deviations themselves. We may call 
it the “identification reach” of the considered model. Accordingly, it should be important 
to try to achieve, for the various models, the maximum sensitivity (i.e., discovery reach) 
as well as highest possible identification reach. 

Here, we will try to quantitatively discuss the above issues in the cases of the Kaluza- 
Klein (KK) graviton exchange in the context of gravity propagating in “large” compactihed 
extra dimensions [DEI, and of the four-fermion contact interactions inspired by the context 
of leptons and quark compositeness HHHj- 1 In particular, our aim will be to assess the 
potential of electron and positron longitudinal polarization at the ILC, to enhance the 
identification reaches on these NP effective interactions. 

Specifically, to this purpose we will take as basic observables the polarized angular 
differential cross sections of the purely leptonic processes, namely, the Bhabha scattering: 

e + + e~ —■> e + + e _ , (1) 

1 Actually, this kind of description more generally applies to a variety of new interactions generated by 
exchanges of very massive virtual objects such as, for example, Z's, leptoquarks and heavy scalars. 
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the annihilation into lepton pairs (l — fi, t): 


e + + e~^l + + r, ( 2 ) 

and the Mpller scattering: 

e~ + e~ —> eT + e~. (3) 

Processes o-o can all receive contributions from graviton exchange and from four- 
fermion contact interactions as well, and represent sensitive probes of the above mentioned 
scenarios, particularly in the case of polarized initial beams [ZHU]. Indeed, the interesting 
feature of the ILC ra is that polarization of both initial beams can be available mi, and 
that in principle this collider could be run in both required modes, nemcly, e + e _ and e _ e _ . 
In Ref. [12], the identification reach on the different contact interactions was studied by 
using a Monte Carlo-based analysis of the unpolarized differential cross sections of processes 
m and 0- Electron and positron longitudinal polarization can substantially help in 
reducing the “confusion” regions of the parameters where models cannot be discriminated 
from each other on a y 2 basis. Accordingly, the identification reach on the individual 
models should be enhanced. 

In Sec. 2 we give the relevant expressions for the differential polarized angular dis¬ 
tributions and for the deviations from the SM predictions in the different NP scenarios 
considered heres. Secs. 3 and 4 outline the numerical analysis and present the numerical 
results for the discovery reaches and the distinction among the models, respectively, as 
allowed by initial beams polarization. Finally, some conclusive remarks are given in Sec. 
5. 


2 Angular distributions and deviations from the SM 

The polarized differential cross section of process m can be conveniently written as (see 
for example [inHE] and P): 

da(P~,P + ) = (1 + P~)(1-P + ) da R | (1-P~)(1 + P+) dcr L 

d cos 9 4 d cos 9 4 d cos 9 

(1 + P-)(1 + P+) da RL , (1 - P~) (1 - P+) da LR , t 

4 d cos 9 ' 4 d cos 9' 1 1 


with 9 the angle between incoming and outgoing electrons in the c.m. frame and P T the 
longitudinal polarization of electron and positron beams, respectively. In Eq. ©: 


where 


d(TLL 

dcos# 


dcr L _ do~ L L da L R,s 

d cos 6 d cos 6 d cos 6 ’ 
dcqt _ da RR da R L ; s 

d cos 9 d cos 9 d cos 9 ’ 


27ra 2 


Gll,s + G 


LL ,t\ 


dcr RR _ 27TQ: e m | ^ ^ \ 2 

3-3 —-yrRR,s + <~T RRt , 

d cos 9 s 1 1 


( 5 ) 
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dg L R,t 

dcos# 


dcr RL ,t _ 27T«g m , ,2 d(JLR, s _ d(jRL,s _ 27TO: 2 m , ,2 

i /j ^LR^ ? i r\ i /j ^LRjS 

d cos d s' d cos 6 d cos 6 s 1 1 


( 6 ) 


The hclicity amplitudes G a p (a,/3 = L,R) can be expressed in terms of the SM 7 and Z 
exchanges in the s- and t-channels, plus deviations due to contact interactions representing 
the physics beyond the SM, as follows: 

° LL ' = “ {~s + + All '*) ’ Gll '‘ = “ (} + rvhf + Ali -‘) ’ 

Gr Rj = u (i + + A RR ,,) , G RR ,« = U (i + + Ariu) , 

+ A L R,4^ • (7) 

Here: u,t = —s(l±cos 6>)/2; M z represents the mass of the Z\ g R = tan 9 W , = — cot 2 6 *^/ 

are the SM right- and left-handed electron couplings to the Z , with 9\v the electroweak 
mixing angle. 

The advantage of Eqs. 0-0 is that the expression for the differential cross section of 
process 0 can be obtained directly from crossing symmetry: one has to replace s <-> w, 
(P~,P + ) —> (Pi , — P 2 ) with Pi and P 2 - denoting the polarizations of the two initial 
electrons, and divide the cross section by 2 to account for identical particles. Also, the 
same equations are easy to adapt to the cross section of the annihilation process 0 , one 
simply must drop all f-channel poles (actually, for this process, in principle 6 lr 7 ^ e RL [THTj). 

Turning to the amplitudes deviations A a p in Eq. 0, in the ADD graviton exchange 
scenario only gravity can propagate in two or more extra spatial dimensions of the 
millimeter order, whereas the SM particles must live only in the ordinary four-dimensional 
space . 2 Massless graviton exchange in extra dimensions translates to the exchange of a 
tower of evenly spaced KK massive states (with vertices given in |17l IH] , see also H3), 
and this effect can be parameterized by the effective, dimension- 8 , contact interaction 
Lagrangian JTj 

£ add = i 4A ( 8 ) 

a h 

Here, is the energy-momentum tensor of the SM particles, Ah is a ultraviolet cut¬ 
off on the summation over the KK spectrum, expected in the (multi) TeV region, and 
A = ±1 (models with A = 1 or A = —1 are denoted as ADD+ and ADD—, respectively). 
The explicit expressions of the corresponding corrections to the SM amplitudes relevant to 
Bliabha scattering are reported in Table UJ In this table, we include also the deviations 
corresponding to models with TeV _1 -scale extra dimensions, parameterized by the “com- 
pactification scale” Me , where also the SM gauge bosons may propagate in the additional 
dimensions [ 2 U E2j. Current experimental limits, from LEP2 and Tevatron, are in the 
range A h > 1.1 — 1.3 TeV J22H2E]- For the TeV -1 -scale extra dimension scenario the limit, 
mostly determined by LEP data, is M c > 6.8 TeV [SJJ. 

2 Actually, while the case of two extra dimensions may be marginal due to cosmological arguments and 
direct gravity experiments, three extra dimensions have been advocated from dark-matter observations 
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Tabic 1: Parametrization of the A Qj g functions in different models for e + e 


new physics model 

^ af3 

composite fermions !5( 

A A 1 

^0.(3,s — ^-ap,t ~ a 2 

a e.m. A a g 

TeV _1 -scale extra dim. [21] [22] 

f 71 " 1 

^a/3,s — ~ (QeQf ± 9a 9g) <1, M 2 

ADD model 12- EH] 

7 7 * 7 

All, s — A RRjS — ,4 i u + a 5 ) 

TTfte.m 4 

. . A . 3 . 

All,* — A rr t — A 4 i u + a t) 

na e . m A% 4 
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A L r, s - .4 (*+ /) 
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The four-fermion contact interaction scenario can be represented by the following 
vector-vector, dimension- 6 , effective Lagrangian (with r/ Q/ g = ±1,0 and a, (5 = L, R) [Sj: 


£ CI 


47T 

1 ± s ef 


E 


T? (e a ^e a ) (//37 m //3) , 
IV a/3 


(9) 


where A a g denote compositeness scales and <5 e / =1 (0) for / = e (/ ^ e). The most popular 
four-fermion contact interaction models (Cl) are defined by specializing, in Eq. ©, the 
helicities according to Table [21 


Table 2: Definition of four-fermion Cl models 


CI model 

Vll 

Vrr 

Vlr 

Vrl 

LL 

±1 

0 

0 

0 

RR 

0 

±1 

0 

0 

LR 

0 

0 

±1 

0 

RL 

0 

0 

0 

±1 

VV 

±1 

±1 

±1 

±1 

AA 

±1 

±1 

±1 

±1 


The corresponding deviations A a/ g, that appear in Eq. (f71) . are listed in Table [U Current 
limits on As significantly vary according the process studied and the kind of analysis 
performed there. In general, the lower bounds are of the order of 10 TeV (a detailed 
presentation can be found in the listings of Ref. JS]). 

It should be noted from Tabled that, contrary to the other cases, for the ADD scenario 
the amplitudes deviations are z-dependent (z = cos 9), and consequently add extra terms, 
proportional to z 3 and z 4 , to the SM angular distribution of the annihilation process 02D- 
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However, these contributions to the total cross section are expected to be tiny, because 
the (expected dominant) interference with the SM amplitudes vanishes when integrated 
over the full angular range and there remain only the terms quadratic in (l8l) . suppressed 
by the corresponding very high power of the small parameter 1/A#. On the contrary, the 
interference between graviton exchange and 1-channel SM exchanges for process (CD may 
give non-vanishing contributions, which may favorably combine with the larger statistical 
precisions expected in this channel and increase the discovery reach on A#. 


3 Discovery reach on the contact interaction models 


We here briefly outline the derivation of the expected discovery reaches on the New Physics 
scenarios introduced in the previous section. The basic objects are the relative deviations 
of observables from the SM predictions due to the NP: 


A(0) 


e>(SM + NP) - 0(SM) 

c>(sm) 


( 10 ) 


and, as anticipated, we concentrate on the differential cross section, O = dcr/dcosfh To 
get an illustration of the effects induced by the individual NP models, we show in Figs. 1-3 
the angular behavior of the relative deviations m for the three leptonic processes under 
consideration (with unpolarized beams), for c.m. energy a /s = 0.5 TeV and selected 
values of the relevant mass scale parameters close to their discovery reaches (unpolarized 
cross sections). The superscript “+” on the Cl mass scales A a p denotes the choice rj a/ 3 = 
1 in Table 0 while the notation ADD± corresponds to A = ±1 in Eq. (JHJ). Vertical 
bars represent the statistical uncertainty in each angular bin, for an integrated luminosity 
£ int = 100 ftV 1 . Of course, the comparison of deviations with statistical uncertainties is an 
indicator of the sensitivity of an observable to the individual effective interaction models. 

Basically, a y 2 analysis of the differential cross section of processes m-m can be 
performed by dividing the angular range into bins and introducing the sum over bins: 



where the relative deviations A (O) are defined in Eq. (FTTH) and SO denotes the expected 
experimental relative uncertainties, that combine statistical and systematic ones. 

As a criterion to constrain the individual models, in particular to set the discovery reach 
on the relevant mass scales, one essentially looks for the smallest values of such parameters 
above which the deviation from the SM prediction is too small to be observable within 
the experimental accuracy. This value, indicated by non-observation of deviations, results 
from the condition 

A 2 < XL, (12) 

where Xcl represents a critical value and we will take x'cl =3.84 for 95% C.L. 

To make contact to the foreseeable experimental situation, we impose cuts in the for¬ 
ward and backward directions of processes Specifically, for Bliabha and Mpllcr 

scattering we consider the cut angular range | cosd| < 0.9 and divide it into nine equal-size 
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-0.8 -0.4 0 0.4 0.8 


cos(0) 

Figure 1: Relative deviations of the unpolarized Bhabha differential cross section from the 
SM prediction as a function of cos 9 at a/s = 0.5 TeV for the Cl models of Table O AA 
(A aa =48 TeV), VV (A+ v =76 TeV), LL (A+ L =37 TeV), RR (A+ R =36 TeV), LR (A+ R =60 
TeV); for the TeV -1 model (M c =12 TeV) and the ADD± models (A^=4 TeV). The 
vertical bars represent the statistical uncertainty in each bin for £ int = 100 ftV 1 . 


bins of width Acos0 = 0.2. Similarly, for the annihilation into muon and tau pairs the 
angular range |cosd| < 0.98 will be considered. To assess the statistical uncertainty, we 
assume the reconstruction efficiency e = 95% for final l + P events (/ = /i, r), and e ~ 100% 
for final e + e _ and e _ e _ pairs. Also, to assess the dependence of discovery reaches on the 
c.m. energy and on the time-integrated luminosity, in the sequel an 1LC with i/s = 0.5 
TeV and 1 TeV will be considered, with £i nt (e + e“) ranging from 100 fb _1 up to 1000 fb -1 . 

Concerning systematic uncertainties, an important source is represented by the uncer¬ 
tainty on beams polarizations, for which we assume 5P ± /P ± = 0.2% (and 5Pb 2 /Pf 2 = 
0.2% for the case of Mpllcr scattering). Also, a systematic error of 0.5% in the luminos¬ 
ity determination is assumed. In the case of the processes 0 and 0 , we analyze the 
following combinations of beams polarizations: (P”,P + ) = (|P~|, —|P + |); (—|P~|, |P + |; 
(|P~|, |P + |); (— |P~|, —|P + |) with the “standard” envisaged values \P~\ = 0.8 and |P + | = 
0.6 (|Pf | = |P 2 - | = 0.8 for Mpller scattering). For the annihilation process 0 , we limit to 
the (P~, P + ) = (|P”|, —|P + |) and (—|P _ |, |P + |) configurations. As for the time-integrated 
luminosity, for simplicity we assume it to be equally distributed between the different po¬ 
larization configurations defined above, and £mt(e _ e _ ) ~ \C mt {e + e~) to account for the 
reduction in luminosity of the e - e - mode due to anti-pinching in the interaction region 
[22' • The discovery limits are derived by taking the sum of the y 2 relevant to the individual 
configurations of polarizations mentioned above and imposing the constraint m- Also, 
we take into account correlations between the different polarized cross sections, but not 
those among individual angular bins. 


7 














0.04 



-0.04 - 1 -*- 1 -*- 1 -*- 1 -*- 1 - 

-0.8 -0.4 0 0.4 0.8 

cos(0) 

Figure 2: Same as in Fig. 1 but for e + e _ —> /i + /r _ , for the Cf models of Table [21 AA 
(A+ a =80 TeV), VV (A+ v =90 TeV), LL (A+ L =45 TeV), RR (A+ R =42 TeV), LR (A+ R =41 
TeV), RL (A rl =43 TeV); for the TeV^ 1 model (M c — 17 TeV) and the ADD± models 
(A#=2.8 TeV). 


Regarding theoretical inputs, for the SM amplitudes we use the effective Born approx¬ 
imation ra taking into account electroweak corrections to propagators and vertices, with 
m top = 175 GeV and mn = 120 GeV. Among the O(a) QED corrections, the numerically 
most important ones come from initial-state radiation. In the case of processes (HD and ED, 
we account for that effect by a using a structure function approach including both hard and 
soft photon emission PH and the flux factor method m, respectively. To minimize the 
effect of radiative flux return to the s-channel Z-exchange, a cut is applied on the radiated 
photon energy, A = < 1 — M|/s with A = 0.9, in order that interactions occur 

close to the nominal collider energy and thus the best sensitivity to the manifestations of 
non-standard physics can be obtained. Other QED effects, such as final-state and initial- 
final state emission, are found in process 0 to be numerically unimportant for the chosen 
kinematical cuts, by using the ZFITTER code • Regarding process ED, the lowest-order 
corrections to the polarized cross sections are evaluated by means of the FORTRAN code 
MOLLERAD pODESI, adapted to the present discussion. 

The numerical results for the discovery reaches on the effective contact interaction 
mass scales, obtained by the above procedure in the case of an ILC with yfs = 0.5 TeV 
and Ant(e + e~) = 100 fb _1 , are summarized for the different processes and polarization 
configurations in Table 01 In this table, only the results for positive interference between 
SM and non-standard interactions are reported (i.e.,r) a p = +1 for Cl models and A = 1 for 
ADD), because the sensitivity reach for negative interference is practically the same. The 
results displayed in Table 01 relevant to the processes CD and ED with unpolarized beams 
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Figure 3: Same as in Fig. 1 but for M 0 ller scattering for the Cl models of Table [21 AA 
(A+ a= 55 TeV), VV (A+ v =55 TeV), LL (A+ L =44 TeV), RR (A+ R =44 TeV), LR (A+ R =21 
TeV); TeV -1 model (Me— 12 TeV) and ADD± models (A^=4 TeV). Here, £ int (e - e - ) ~ 
|£i nt (e + e _ ) has been assumed. 


are consistent with the recent estimates in Ref. ED- Also, the limits on the Cl mass scales 
resulting from l + l~ final states are based on the assumption of fi — t universality. 

Some features of the numbers in Table |21 are noteworthy. Firstly, there is complemen¬ 
tarity of the leptonic processes under consideration, ©■ © and in the search for Cl 
and ADD scenarios. Secondly, for the chosen values of the time-integrated luminosity, the 
discovery reaches on A h of the Bhabha and Mpller scattering processes can be larger than 
~ 8 ^. Also, Bhabha scattering is the process leading to the best search reach for A#, 
due to the available higher statistics in the data sample in comparison to either Mpller 
scattering or, to a large extent, to the muon pair production process. 

One may recall that, on the purely statistical basis, the discovery reach on A h should 
be expected to scale with s and £j nt like ~ (s'^-Cint) 1 / 8 , and the discovery reach on the Cl 
mass scales As like ~ (s/lint) 1 / 4 , according to the different dimensions of the operators © 
and ©. 


4 Distinction among the New Physics models 

Let us assume one of the models to be consistent with data and call it “true” model, 
for example the ADD model © with some value of A h- We want to assess the level at 
which this “true” model is distinguishable from the other ones, that can compete with 
it as sources of corrections to the SM and we call “tested” models, for any values of the 
corresponding mass scale parameters. For example, we may take as “tested” model any 
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Table 3: Discovery reach (in TeV) on the mass scale parameters (95% C.L.) from the 
lepton pair production processes at y/s = 0.5 TeV. For the e + e _ mode the three entries 
refer to £i nt (e + e _ ) = 100 / 6 _1 and the polarizations configurations (|P _ |, |P + |) = ( 0 , 0 ); 
(0.8,0); (0.8,0.6). For the e~e~ mode the configurations are (|Pf|, |P 2 _ |) = (0,0); (0.8,0); 
( 0 . 8 , 0 . 8 ) and £ int (e _ e _ ) « ~£ int (e + e~). 


model 

process 

e + e 

" — > e 

+ e~ 

e e~ 

" —*■ e 

e 

e + e 

-»• £ 

+ h~ 

e + e 

- -> i+r 

ADD± 

(A h) 

4.1 

4.2; 

4.3 

3.8 

4.0; 

4.1 

2.8 

2 . 8 ; 

2.9 

3.0 

3.0; 

3.2 

VV 

(A) 

76.2; 

80.8; 

86.4 

64.0; 

68 . 8 ; 

71.5 

75.5; 

76.4; 

83.7 

OO 

QO 

90.7 

99.4 

AA 

(A) 

47.4; 

49.1; 

69.1 

58.0; 

62.0; 

64.9 

67.3; 

68 . 2 ; 

74.8 

80.1; 

81.1 

88.9 

LL 

(A) 

37.3; 

45.5; 

52.5 

43.9; 

52.4; 

55.2 

45.0; 

51.0; 

57.5 

53.4; 

60.5 

68.3 

RR 

(A) 

36.0; 

44.7; 

52.2 

42.3; 

52.3; 

55.4 

43.2; 

50.6; 

57.5 

51.3; 

60.0 

68.3 

LR 

(A) 

59.3; 

61.6; 

69.1 

20 . 1 ; 

22 . 1 ; 

31.5 

40.6; 

46.0; 

52.6 

48.5; 

55.0 

62.8 

RL 

(A) 

A RL — A LR 

Arl — A 

LR 

40.8; 

46.7; 

53.4 

48.7; 

55.6 

63.6 

TeV 

(M c ) 

12 . 0 ; 

12 . 8 ; 

13.8 

11.7; 

12.5; 

12.9 

16.8; 

17.1; 

18.7 

20 . 0 ; 

20.3 

22.2 


one of the four-fermion effective contact interactions listed in Table |2] To that purpose, we 
can introduce relative deviations of the differential cross sections from the ADD predictions 
in each angular bin, arising from the Cl models, analogous to Eq. (HDD: 


A(0) 


e>(ci) - c>(add) 

£>(ADD) 


(13) 


Correspondingly, a y 2 function analogous to Eq. m can be introduced, with 8(0) defined 
in the same way as 8(0) but, in this case, the statistical uncertainty is referred to the ADD 
model and therefore depends on the particular value of A r. Since there can be “confusion” 
regions of A h and A a p values where also some Cl model can be consistent to the ADD 
predictions, on the basis of such y 2 we can study whether these “tested” models can be 
excluded or not to a given confidence level, that we always assume 95%, once the ADD 
model has been assumed as “true”. Then, we scan all values of A H up to the discovery 
reach. 

Thus, let us choose anyone of the “tested” Cl models in m, for definiteness the 
VV one defined in Table 0 so that the y 2 mentioned above will be a function of the 
two parameters eci = r]/A and 6r = A/A# defined in Eqs. © and (JHJ), respectively. 
In Fig. El the four gray areas in the (eci, €h) plane, corresponding to the sign choices 
(rj, A) = (1,1); (-1,-1); (1,-1); (—1,1), represent values of the parameters for which 
both the ADD and the VV models can give observable effects in unpolarized Bhabha 
scattering at the ILC, with 95% C.L. Also, the horizontal and vertical bands correspond to 
the discovery reaches on the ADD and VV models at the 95% C.L., derived in the previous 
section in the unpolarized case. The “confusion region” is the area where the y 2 is smaller 
than Xcl =3.84 and the two models cannot be distinguished at the 95% C.L. 

As indicated in Fig. El one can find a maximal absolute value of the scale parameter e# 
for which the “tested” VV model hypothesis is expected to be excluded at the 95% C.L. 
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Figure 4: Region of confusion at 95% C.L. for ADD and VV models obtained from unpo¬ 
larized Bhabha scattering at yfs =0.5 TeV and £ in t(e + e _ ) =100 fb -1 . 


for any value of the Cl parameter eci taken in the gray area. We denote the corresponding 
ADD mass scale parameter as A)) v and call it “exclusion reach” of the VV model. 

It is worth noticing that the “confusion region” is located only in the areas with ( 77 , A) = 
(1,-1) and (—1,1). This is so because in these regions the ADD± and VV=F models 
provide the same sign of the leading interference with the SM, hence of the deviation 
in the differential cross sections as depicted in Fig. 1 and, therefore, they can mimic each 
other. Conversely, there is no confusion region in the areas with (17, A) = (1,1) and (—1, —1) 
where interferences in the angular distribution have opposite signs and therefore can easily 
be distinguished. 

Also, one may point out that, although this kind of y 2 analysis could in principle be 
applied to any observable, the choice of the differential distribution as basic observable is 
rather crucial to identify the ADD scenario. Integrated observables such as, for example, 
the total cross section, the forward-backward asymmetry and the left-right asymmetry, do 
not allow to derive such a compact confusion region concentrated around the origin as 
in Fig. 01 For instance, the confusion regions determined from the integrated observables 
would extend up to a few units of A/A# (TeV -4 ) and the corresponding reaches on A h 
would be too much below the current discovery limits. The procedure outlined above can 
be repeated for all other types of effective contact interaction models in Table [21 as well 
as the TeV -1 gravity model mentioned in Table d and consequently one can evaluate 
the corresponding “exclusion reaches” A^ A , A^ R , A^ L , A^ R and Ajf v . The results of this 
kind of analysis for the three processes of interest here with unpolarized beams (gray 
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Figure 5: Exclusion and identification reaches on A# at 95% C.L. obtained from the 
leptonic processes at y/s =0.5 TeV, £ int (e + e _ ) =100 fb _1 for unpolarized beams (gray 
histograms) and ~ |£i nt (e + e _ ) and both beams polarized (black histograms) 

are illustrated. 


histograms) as well as polarized beams (black histograms), are represented in Fig. 0 In 
this figure, we have considered an ILC with = 0.5 TeV, time-integrated luminosity 
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Table 4: 95% C.L. identification reach on A® at a/s = 0.5 TeV and £; n t(e + e ) = 100 fb 1 . 


process 

p- = P+ = 0 
A! D (TeV) 

p-| = 0.8; l-P+l = 0.6 
AS? (TeV) 

e + e —»e + e 

2.5 

3.1 

e + e- —> /i + yU 

2.5 

2.6 

e + e“ —> /+/- 

2.7 

2.8 

^int(^ ^ ~ 3 ^int(^ ^ ) 

e~e~ —> e~e~ 

Pi 2 = 0 

o 

bo 

1.4 

2.1 


^int(e + e _ ) = 100 fb- 1 and, in the polarized case, the values of the beams polarizations 
anticipated in the previous section, i.e.: |P _ | = 0.8 and |P + | = 0.6 (and for process © 
one third of the e + e~ luminosity and Pfi = |P 2 _ | = 0.8). 

As the final step, the “identification reach” on the ADD scenario can be defined as the 
minimum of the A H “exclusion reaches”, A® = min{ A)) v , A^ a ,A^ r , A^ l , A^ r , A]) cV } as 
indicated in Fig.0 It is clear that taking A h < A® allows one to exclude all composite-like 
Cl models as well as the TeV _1 -scale gravity model. In the specific example of the ILC 
parameters assumed in the derivation of the results of Fig. El it turns out that the “identifi¬ 
cation reach” on the ADD graviton exchange model obtained from unpolarized (polarized) 
Bliabha scattering corresponds to A® =2.5 (3.1) TeV. The identification reaches on A h 
obtained from the different leptonic processes are summarized in Tabic 01 

The role of beam polarization in increasing the sensitivity of the leptonic processes to 
the graviton exchange effects, and particularly in enhancing the potential of distinction 
from the other NP scenarios, is essential, as exemplified in Fig. El in the case of process (HD- 
Indeed, this figure indicates that the reduction of the “confusion region” with respect to 
the unpolarized case, allowed by beams longitudinal polarization, can be quite substantial. 

Fig. □ shows the 95% C.L. identification reach on the graviton exchange mass scale 
An as a function of time-integrated luminosity 3 and two values of the c.m. energy,= 
0.5 TeV and 1 TeV. The numerical procedure to obtain those results is the y 2 analysis 
described above, separately applied to the individual processes ©%©• In this figure, 
as well as in the subsequent ones, curves are labelled by the final states of the different 
processes, i.e., by e + e - for Bhabha scattering, l + l~ for the combination of and re¬ 

production in the annihilation process, and e _ e- for Mpller scattering. The results for 
the unpolarized differential cross sections are also shown for comparison. In the case of 
polarization, the polarized cross sections have been combined similar to the procedure 
followed in the derivation of the discovery reaches, outlined in the previous section and 
summarized in Table El with the values of longitudinal polarizations exposed there. 

The identification reach on the TeV _1 -scale gravity model, obtained by applying the 
same kind of analysis to the three lepton production processes of interest here, is presented 
in Fig. El while the identification reaches on the four-fermion contact interactions © are 

3 On the horizontal axis the luminosity in the e + e“ channel is reported. One should recall that 1/3 of 
that luminosity is assumed for the e~e~ mode. 
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Figure 6: Region of confusion for ADD and VV models obtained from unpolarized (dashed 
curve) and polarized (solid curve) Bhabha scattering at y/s =0.5 TeV and £ int =100 fb _1 . 

shown in Figs. IH1 and ITU1 In these figures, curves corresponding to values of identification 
reaches that are found to fall below the currently available experimental limits (listed in 
sec. El) have not been included. 

5 Concluding remarks 

We here briefly comment on the main features of the findings in FigsfTHTTH for the identi¬ 
fication reaches on the effective non-standard interactions considered in Sec. El obtainable 
from lepton pair production at the ILC with longitudinally polarized beams. 

Fig. □ shows that, of the three considered processes, Bhabha scattering definitely has 
the best identification sensitivity on the scale Ah characterizing the ADD model of Eq. (JHJ) 
for gravity in “large” compactihed extra dimensions. As one can see, in the polarized 
case, the identification reach ranges from 3.1 TeV to 6.9 TeV, depending on c.m energy 
and on luminosity. It could be of some interest in this regard to estimate the “resolving 
power” on A h obtainable from the polarized process 0. This is defined as the precision 
that in principle might be achieved on the determination of A h, in the case where the 
effects of KK graviton exchange were observed. Fig. im shows the uncertainty obtainable 
on A h when we vary this parameter in the range between the current experimental bound 
and the expected identification reaches. Of course, the larger A h the worse the precision 
on it. In contrast, although being competitive with the other processes for the discovery 
reach, Mpller scattering is found to give very poor identification reaches on A h even in 
the polarized case, and appears to be saturated by the systematic uncertainty as indicated 


e e e e 
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Figure 7: 95% CL identification reach on the cutoff scale A h in the ADD model as a 
function of the integrated luminosity obtained from lepton pair production processes with 
unpolarized and both polarized beams at ILC(0.5 TeV) (left panel) and ILC(1 TeV) (right 
panel). 


by the corresponding curves quickly becoming flat for increasing luminosity. Indeed, for 
this for process the KK graviton exchange contribution to the RR and LL cross sections in 
Eq. © is cos ^-independent therefore identical to the Cl interactions ©, while the cos 9- 
dependent contribution to the LR cross section is considerably suppressed kinematically. 
Accordingly, the potential of identification from the angular distributions is reduced. This 
is qualitatively true in general, i.e., the identification reaches from process © are rather 
moderate, and for some models well-below the discovery reach. 



Figure 8: Same as in Figd but for the compactihcation scale Me- 
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Figure 9: Same as in FigjTJbut for the compositeness scale for VV (top panel), AA (central 
panel) and LL (lower panel) type interactions. 
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Figure 10: Same as in Fig[7]but for the compositeness scale for RR (top panel), LR (central 
panel) and RL (lower panel) type interactions. 
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Figure 11: Resolving power 5Ah on the mass scales (la level) vs. A h obtained from the 
polarized Bhabha scattering for the c.m. energy a/s =0.5 TeV (left panel) and 1 TeV (right 
panel) and integrated luminosity of £ int = 100 ftV 1 and 1000 fb _1 . 


The l + l~ annihilation process (J2J) definitely provides the most stringent identification 
reaches on TeV _1 -scale gravity in extra dimensions and on the four-fermion contact inter¬ 
actions m as well. Quantitatively, Figs. IHUTTil show that the identification reach on the 
compactihcation mass parameter Me ranges from 15 TeV to 35 TeV depending on the ILC 
energy and luminosity, and that the reaches on Cl four-fermion interactions © are rather 
high and significantly enhanced by polarization: roughly, depending on the energy and the 
luminosity, the identification reach on Ayy ranges between 62 TeV and 160 TeV; that on 
A A a between 70 TeV and 170 TeV; for A LL between 55 TeV and 135 TeV. Finally, for A RR , 
Al R and A R l the discrimination range can be between 57 TeV and 142 TeV. 

As regards the effective role of polarization in the derivation of the identification reaches 
in Figs. mm in particular that of the positron beam polarization in addition to the electron 
one, it turns out numerically that unpolarized beams do not lead to sensible “identification 
reaches” on the mass scales relevant to the LL and RR four-fermion Cl models in the case 
of Bhabha and Mpllcr scattering, and on the mass scales of LL, RR, LR, RL models in the 
case of annihilation into l + l~. This relates to the fact that |g£| ~ |g R | so that the (leading) 
interferences of those NP interactions with the SM have almost identical, hence in practice 
indistinguishable, angular dependence (see Eq. ©). Also, unpolarized Mpller scattering 
may give identification limits at the level of current experimental limits, for selected values 
of energy and luminosity. In general then, for a given process, polarized beams can have 
much higher identification potential than the unpolarized ones. Specifically, it turns out 
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that the identification reaches in Figs. [7|EI]for LL and RR models from Bhabha scattering 
and those on LL, RR, LR, RL and gravity in TeV -1 extra dimensions from annihilation into 
l + l~ , are mostly rely on electron polarization and not as much on the simultaneous positron 
polarization. By contrast, the identification reaches on ADD, AA, VV, LR and TeV -1 
models from Bhabha scattering and those on ADD, AA and VV models from annihilation 
into l + l~ pairs definitely need both electron and positron polarization. It may also be 
noticed that, in many cases, polarized cross sections allow to identify the Cl models almost 
near the discovery reaches (summarized for the lowest considered values of energy and 
luminosity in Table 0 ). 

Acknowledgements 

AAP acknowledges the support of INFN and of MIUR (Italian Ministry of University and 
research). NP has been partially supported by funds of MIUR and of the University of 
Trieste. This research has been partially supported by the Abdus Salam ICTP and the 
World Federation of Scientists (National scholarship programme). AAP and AVT also 
acknowledge the Belarusian Republican Foundation for Fundamental Research. 


References 

[1] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) 

| arXiv:hep-ph/98033151. 

[2] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Rev. D 59, 086004 (1999) 

| arXiv:hep-ph/9807344]]. 

[3] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 
257 (1998) | arXiv:hep-ph/98043981. 

[4] G. ’t Hooft, in “Recent Developments In Gauge Theories”, Proceedings, Nato Ad¬ 
vanced Study Institute, Cargese, France, August 26 - September 8, 1979, edited by 
G. ’t Hooft, C. Itzykson, A. Jaffe, H. Lehmann, P. K. Mitter, I. M. Singer and R. Stora 
New York, USA: Plenum (1980) (Nato Advanced Study Institutes Series: Series B, 
Physics, 59)] 

S. Dimopoulos, S. Raby and L. Susskind, Nucl. Phys. B 173, 208 (1980). 

[5] E. Eichten, K. D. Lane and M. E. Peskin, Phys. Rev. Lett. 50, 811 (1983). 

[6] R. Riickl, Phys. Lett. B 129, 363 (1983). 

[7] A. A. Pankov and N. Paver, Phys. Rev. D 72, 035012 (2005) JarXiv : liep-ph / 050 11 70 . 

[8] P. Osland, A. A. Pankov and N. Paver, Phys. Rev. D 68 , 015007 (2003) 
j arXiv:hep-ph /03041231. 

[9] A. A. Pankov and N. Paver, Eur. Phys. J. C 29, 313 (2003) arXiv:hep-ph/0209058 . 


19 





[10] J. A. Aguilar-Saavedra et al. [ECFA/DESY LC Physics Working Group Collabora¬ 
tion], “TESLA Technical Design Report Part III: Physics at an e + e“ Linear Collider,” 
DESY-01-011, arXiv:hep-ph/0106315; 

T. Abe et al. [American Linear Collider Working Group Collaboration], “Linear 
collider physics resource book for Snowmass 2001. 1: Introduction,” in Proc. of the 
APS/DPF/DPB Summer Study on the Future of Particle Physics (Snowmass 2001) 
ed. N. Graf, SLAC-R-570, arXiv:hep-ex/0106055 

K. Abe et al. [AGFA Linear Collider Working Group Collaboration], 
arXiv:hep-ph/0109166 

[11] G. Moortgat-Pick et al., arXiv:hep-ph/0507011 

[12] G. Pasztor and M. Perelstein, in Proc. of the APS/DPF/DPB Summer Study on the 
Future of Particle Physics (Snowmass 2001) ed. N. Graf, arXiv:hep-ph/0111471 

[13] B. Schrempp, F. Schrempp, N. Wermes and D. Zeppenfeld, Nucl. Phys. B 296, 1 
(1988). 

[14] F. Cuypers, eConf C960625, NEW137 (1996) arXiv:hep-ph/9611336|. 

[15] S. Cullen, M. Perelstein and M. E. Peskin, Phys. Rev. D 62, 055012 (2000) 
|arXiv:hep-ph/0001166|. 

[16] B. Qin, U. L. Pen and J. Silk, arXiv:astro-ph/0508572. 

[17] G. F. Giudice, R. Rattazzi and J. D. Wells, Nucl. Phys. B 544, 3 (1999) 

jarXiv:hep-ph/981129lJ. 

[18] T. Han, J. D. Lykken and R. J. Zhang, Phys. Rev. D 59, 105006 (1999) 

|arXiv:hep-ph/9811350|. 

[19] G. F. Giudice, T. Plehn and A. Strumia, Nucl. Phys. B 706, 455 (2005) 

| arXiv:hep-ph/04083201. 

[20] J. L. Hewett, Phys. Rev. Lett. 82, 4765 (1999) arXiv:hep-ph/9811356|. 

[21] K. M. Cheung and G. Landsberg, Phys. Rev. D 65, 076003 (2002) 

| jarXiv:hep-ph /01103461. 

[22] T. G. Rizzo and J. D. Wells, Phys. Rev. D 61, 016007 (2000) |arXiv:hep-ph/9906234 . 

[23] S. Ask, arXiv:hep-ex/0410004. 

[24] M. K. Unel [for the CDF and DO Collaborations], FERMILAB-CONF-04-206-E, 
arXiv: hep-ex/0411067 

[25] G. Landsberg [CDF Collaboration], eConf C040802, MOT006 (2004), 

(arXiv: hep-ex/0412028]|. 

[26] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 95, 161602 (2005) 

| arXiv:hep-ex/05060631. 


20 








[27] For a review see, e.g., K. Cheung, arXiv:hep-ph/0409028, 

[28] S. Eidelman et al. [Particle Data Group], Phys. Lett. B 502, 1 (2004). 

[29] J. E. Spencer, lilt. J. Mod. Phys. A 11, 1675 (1996). 

[30] M. Consoli, W. Hollik and F. Jegerlehner, CERN-TH-5527-89 Presented at Workshop 
on Z Physics at LEP ; 

G. Altarelli, R. Casalbuoni, D. Dominici, F. Feruglio and R. Gatto, Nucl. Phys. B 
342, 15 (1990). 

[31] For reviews see, e.g., O. Nicrosini and L. Trentadue, in Radiative Corrections for 
e + e~ Collisions, ed. J. H. Kiihn 25 (Springer, Berlin, 1989), p. 25; in QED Structure 
Functions, Ann Arbor, MI, 1989 , ed. G. Bonvicini, AIP Conf. Proc. No. 201 (AIP, 
New York, 1990), p. 12. 

[32] For a review see, e.g., W. Benakker and F. A Berends: Proc. of the Workshop on 
Physics at LEP2 , CERN 96-01, vol. 1, p. 79 and references therein. 

[33] D. Bardin, P. Christova, M. Jack, L. Kalinovskaya, A. Olchevski, S. Ricmann and 
T. Riemann, Comput. Phys. Commun. 133, 229 (2001) hep-ph/9908433|. 

[34] V. A. Mosolov, N. M. Shumeiko and J. G. Suarez, Int. J. Mod. Phys. A 15 (2000) 
2377. 

[35] N. M. Shumeiko and J. G. Suarez, J. Phys. G 26, 113 (2000) |arXiv:hep-ph/9912228 . 

[36] A. Ilyichev and V. Zykunov, Phys. Rev. D 72, 033018 (2005) |arXiv:hep-ph/0504191 . 

[37] D. Bourilkov, arXiv:hep-ph/0305125 


21 




